Device for calculating positional data of standard points of photogrammetric target

ABSTRACT

A three points range of high brightness and a four points range of high brightness in a photographed image are selected. It is checked if an identical end point is shared by the three and four points range selected. It is checked if a center of gravity corresponding to the identical end point, a center of gravity corresponding to another end point of the three points rage opposite to the identical end point, and a center of gravity corresponding to another end point of the four points range opposite to the identical end point, meet a positional relationship between standard points of a target.

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to a calculation of camera position(photographing position) in a photogrammetric analytical measurementsystem, in which a survey map is produced on the basis of a set ofphotographed pictures obtained at two different photographing positions.

2. Description of the Related Art

For example, photogrammetry is carried out at a traffic accident spot.The traffic accident spot is photographed by a camera in at least twodifferent positions, and a survey map of the traffic accident spot isproduced based on a set of photographed pictures obtained at thedifferent positions.

In particular, a two-dimensional coordinate system is defined on each ofthe photographed pictures, and two-dimensional positions of the objects,which are recorded on each picture, are determined by thetwo-dimensional coordinate system. Then, a three-dimensional coordinatesystem is defined on the basis of the two sets of two-dimensionalcoordinate systems, and three-dimensional positions of the recordedobjects are determined from the three-dimensional coordinate system.Accordingly, it is possible to produce a survey map of the trafficaccident spot by drawing the objects on a sheet of paper in such amanner that the objects are projected on one of the three planes definedby the three-dimensional system.

Before accurately scaled distances and lengths can be reproduced on thesurvey map, a standard measurement scale must be recorded together withthe objects in the photographed pictures. Also, a standard referenceplane, on which the survey map should be drawn, must be defined in thephotographed pictures.

Usually, in order to define the standard measurement scale and thereference plane, three respective cone-shaped markers, which areidentical to each other, are positioned at suitable locations around atraffic accident spot. Namely, a distance between two apexes of thecone-shaped markers is measured, for example, with a measuring tape, andset as the standard measurement scale, and a plane, defined by the threeapexes of the cone-shaped markers, is utilized as the reference plane.As the positioning of the cone-shaped markers and the measurement of thedistance between the two apexes are performed by operators, a process ofpreparing for photographing is onerous and requires a lot of time.

In order to solve the above mentioned problems, photogrammetric analyticmeasurement, using a target which has a triangular frame instead of thethree cone-shaped markers, is disclosed in Japanese Unexamined PatentPublication (Kokai) No. P10-221072. Standard point members arerespectively mounted on each of three apexes of the target. A distancebetween the standard point members is set as the standard measurementscale, and a plane defined by the standard point members is utilized asthe reference plane. A reflecting sheet, for example, is attached on thestandard point members so that the standard point members can be viewedeasily in photographed pictures.

The determination of the three-dimensional positions of the recordedobjects from the three-dimensional coordinate system is carried out byiterating a series of approximate calculations, using a computer havinga monitor on which the set of photographed pictures is displayed.

Initially, before a first series of approximate calculations isexecuted, the standard point members and a suitable point on one of therecorded objects in the set of photographed pictures displayed on themonitor, are selected with a cursor, by manipulating a mouse, wherebythree sets of two-dimensional coordinates, corresponding to the standardpoint members, and a set of two-dimensional coordinates on the suitablepoint of the selected object are inputted to the computer. Namely, theexecution of the first series of approximate calculations is based uponthe inputted two-dimensional coordinates of the standard point membersand the inputted two-dimensional coordinates of the suitable point onthe selected object, thereby determining three-dimensional positions ofthe standard point members and a three-dimensional position of thesuitable point of the selected object from the three-dimensionalcoordinate system.

Then, another suitable point on the selected object in the set ofphotographed pictures displayed on the monitor is indicated with thecursor, by manipulating the mouse, whereby a set of two-dimensionalcoordinates of the other suitable point on the selected object isinputted to the computer, and a second series of approximatecalculations is executed on the basis of the inputted two-dimensionalcoordinates of the other suitable point on the selected object, therebydetermining a three-dimensional position of the other suitable point onthe selected object from the three-dimensional coordinate system. Thisprocedure is continued until a sufficient number of points on theselected object are indicated, to thereby specify a three-dimensionalprofile of the selected object.

The aforementioned series of approximate calculations is executed basedon positional data of the camera at which the pictures are photographed.The positional data of the camera includes a distance from the camera tothe selected object, an angle of the camera against the object. Thepositional data of the camera is calculated based on the two-dimensionalcoordinates of the standard point members selected by the operator.

However, it takes a comparatively long time to select the standard pointmembers in the set of photographed pictures displayed on the monitor,with the cursor, by manipulating the mouse. Therefore, when there are alot of the photographed pictures, it becomes onerous for the operator tooperate the above mentioned manipulation of the mouse. Further, theaccuracy of the photogrammetric analytic measurement is dependent uponthe skill of the operator.

SUMMARY OF THE INVENTION

Therefore, an object of the present invention is to provide a device forcalculating positional data of standard points of a target ofphotogrammetric analytic measurement.

In accordance with an aspect of the present invention, there is provideda device, for calculating positional data of standard points of a targetof photogrammetric analytic measurement, comprising an image dataobtaining processor and a value of two-dimensional coordinatescalculating processor.

The image data obtaining processor photographs an object ofphotogrammetric analytic measurement with a target that has a firststandard point and a second standard point and a third standard point.

Each length of distance between the first standard point and the secondstandard point and the third standard point is predetermined. An angledefined by a first straight line connecting the first standard point andthe second standard point, and a second straight line connecting thesecond standard point and the third standard point, is predetermined. Anumber of assistant points provided on the first straight line and anumber of assistant points provided on the second straight line aredifferent.

The value of two-dimensional coordinates calculating processorcalculates each two-dimensional coordinate value of the first standardpoint, the second standard point and the third standard point in aphotograph coordinate system that is a two-dimensional coordinate systemof a photographed image, obtained by the image data obtaining processor,after determining each position of the first standard point, the secondstandard point and the third standard point in the photographed image,by comparing the number of the assistant points provided on the firststraight line and the number of the assistant points provided on thesecond straight line.

The first standard point, the second standard point, the third standardpoint and the assistant points provided on the first straight line andprovided on the second line are made of a reflecting material, such thatthe first standard point, the second standard point, the third standardpoint and the assistant points have a relatively higher brightness thanother objects in the photographed image.

The value of two-dimensional coordinates calculating processor isprovided with: a binarization processor that classifies each pixelcomposing of digital image data of the photographed image into a highbrightness class and a low brightness class after comparing a brightnessof each pixel with a predetermined threshold value; a group extractingprocessor that extracts groups in which a plurality of the pixels of thehigh brightness class sequentially exist in the digital image data, anarea of which is within a predetermined range; a calculating processorthat calculates a value of two-dimensional coordinates of a center ofgravity of each of the groups, based on the brightness of each pixelcomposing the groups and a value of coordinates of the each pixel in thephotograph coordinate system; a first points range extracting processorthat extracts a first points range on which a plurality of the centersof gravity lie, a number of the plurality of the centers of gravity onthe first points range being equal to a total point number including thefirst standard point, the second standard point and the number of theassistant points on the first straight line; a second points rangeextracting processor that extracts a second points range on which aplurality of the centers of gravity lie, a number of the plurality ofthe centers of gravity on the second points range being equal to a totalpoint number including the second standard point, the third standardpoint and the number of the assistant points on the second straightline; a combination extracting processor that extracts a combination ofthe first points range and the second points range in which the firstpoints range and the second points range share the center of gravity ofan identical end point; and a standard point determining processor thatdetermines a center of gravity of an end point of the first points rangeselected by the selecting processor, being opposite to the identical endpoint shared by the first points range and the second points range, asthe first standard point, determines the center of gravity of theidentical end point shared by the first points range and the secondpoints range as the second standard point, and determines a center ofgravity of an end point of the second points range selected by theselecting processor, being opposite to the identical end point shared bythe first points range and the second points range, as the thirdstandard point.

The number of the assistant points provided on the first straight lineis two, and the number of the assistant points provided on the secondstraight line is one, so that the first points range is a four pointsrange in which four centers of gravity range on a same straight line inthe photographed image, and the second points range is a three pointsrange in which three centers of gravity range on a same straight line inthe photographed image.

In accordance with another aspect of the present invention, there isprovided a method for calculating positional data of standard points ofa target of photogrammetric analytic measurement comprising: a firststep that photographs an object of photogrammetric analytic measurementwith a target that has a first standard point and a second standardpoint and a third standard point, by a camera; and a second step thatcalculates each two-dimensional coordinate value of the first standardpoint, the second standard point and the third standard point in aphotograph coordinate system that is a two-dimensional coordinate systemof a photographed image photographed by the camera.

Each length of distance between the first standard point and the secondstandard point and the third standard point is predetermined. An angle,defined by a first straight line connecting the first standard point andthe second standard point, and a second straight line connecting thesecond standard point and the third standard point, is predetermined. Anumber of assistant points provided on the first straight line and anumber of assistant points provided on the second straight line aredifferent.

In the second step, each two-dimensional coordinate value of the firststandard point, the second standard point and the third standard pointin the photograph coordinate system, after determining each position ofthe first standard point, the second standard point and the thirdstandard point in the photographed image based on a difference betweenthe number of the assistant points provided on the first straight lineand the number of the assistant points provided on the second straightline.

In accordance with another aspect of the present invention, there is aprovided a recording medium in which a calculating program of standardpoints of a target of photogrammetric analytic measurement is stored.

The target has a first standard point and a second standard point and athird standard point. Each length of distance between the first standardpoint and the second standard point and the third standard point ispredetermined. An angle, defined by a first straight line connecting thefirst standard point and the second standard point, and a secondstraight line connecting the second standard point and the thirdstandard point, is predetermined. A number of assistant points providedon the first straight line and a number of assistant points provided onthe second straight line are different.

The calculating program comprises: a values of two-dimensionalcoordinate calculating routine that calculates each two-dimensionalcoordinate value of the first standard point, the second standard pointand the third standard point in a photograph coordinate system that is atwo-dimensional coordinate system of a photographed image, photographedin such a manner that an object of photogrammetric analytic measurementis photographed with the target by a camera, after determining eachposition of the first standard point, the second standard point and thethird standard point in the photographed image based on a differencebetween the number of the assistant points provided on the firststraight line and the number of the assistant points provided on thesecond straight line.

BRIEF DESCRIPTION OF THE DRAWINGS

The object and other objects of the present invention will be betterunderstood from the following description, with reference to theaccompanying drawings, in which:

FIG. 1 is a plan view of a target in an operational state, to which anembodiment, according to the present invention, is applied;

FIG. 2 is a side view of the target of FIG. 1;

FIG. 3 is a sectional view of the target of FIG. 1;

FIG. 4 is a plan view of a non-reflecting member viewed from a bottomside of a second bar;

FIG. 5 is an enlarged view of a control unit box of the target;

FIG. 6 is a sectional view of the control unit box of FIG. 5;

FIG. 7 is a conceptual view of a positional relationship between thetarget and a camera;

FIG. 8 is a photographed image photographed by the camera;

FIG. 9 is a conceptual view of a positional relationship between thephotographed image and the target;

FIG. 10 is a conceptual view of a relationship between a scenecoordinate system and a camera coordinate system;

FIG. 11 is a flowchart indicating a process of calculating a position ofa camera;

FIG. 12 is a flowchart of procedure of calculating a center of gravityof a group;

FIG. 13 is a conceptual view of each pixel of a digital image data ofthe photographed image;

FIG. 14 is a first half of a flowchart of procedures of extractioncombinations of the groups;

FIG. 15 is shows the latter half of the flowchart of the procedures ofextraction combinations of the groups;

FIG. 16 is a flowchart of procedures of calculation of a cameraposition; and

FIG. 17 is a block diagram of a coordinate-calculator system, in which aphotogrammetric analytic measurement is performed, according to thepresent invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 1 is a plan view of the target 10, with portions broken away forclarity, and FIG. 2 is a side view of the target 10. The target 10 hasan L-shaped figure, comprising a first bar 12 and a second bar 14. Thefirst and second bars 12, 14 are made of metal material. The first andsecond bars 12, 14 respectively have a shape of a quadratic prism, whichis hollow. A non-reflecting sheet, which does not reflect light, isattached to the whole outer side surface of each of the first and secondbars 12, 14. The width of each of the first and second bars 12, 14 is alength L_(W), and the thickness of each of the first and second bars 12,14 is a length L_(H).

An adhesive is applied to the surface of the non-reflecting sheet, whichcontacts the surface of the bars 12, 14. The other surface of the sheetis colored black, and rough. Incident light on the black and roughsurface is absorbed and diffused, so that a luminance amount ofreflected light is extremely reduced. Note that, for example, a blackflattering or matter agent can be applied on the outer surface of eachof the first and second bars 12, 14, instead of utilizing thenon-reflecting sheet.

A controlling unit box 20, which is shaped parallelepiped, is unitarilyfixed at one end of the first bar 12. The controlling unit box 20 ismade of metal material. The non-reflecting sheet is attached to thewhole outer surface of the controlling unit box 20. The thickness of thecontrolling unit box 20 equals the thickness L_(H) of the first bar 12.The width of the controlling unit box 20 is twice the width L_(W) of thefirst bar 12. The controlling unit box 20 is positioned in such a mannerthat a side surface 20 b and a side surface 12 b of the bar 12 lie on asame plane. A side surface 20 c of the controlling unit box 20 isparallel to a side surface 12 c of the bar 12, being offset in adirection opposite to the side surface 20 b.

One end 14 a of the second bar 14 is rotatably mounted on the sidesurface 20 c of the controlling unit box 20, by a hinge 15. A sidesurface 14 b of the second bar 14 and an end surface 20 d of thecontrolling unit box 20, opposite to the end surface at which the firstbar 12 is fixed, lie on a same plane, when the target 10 is utilized forthe aforementioned photogrammetry.

As shown in FIG. 1, the side surface 14 c of the second bar 14 and theside surface 20 c of the controlling unit box 20 define an angle θ.Namely, the angle θ is a right angle made by an axis 12 x (representedas a broken line) of the first bar 12 and an axis 14 x (represented as abroken line) of the second bar 14. A stay 16, which is a fixing member,is connected to the first and second bars 12, 14, at the side of theangle θ. Rotational movement of the second bar 14 is prevented by thestay 16, so that the positional relationship between the first and thesecond bars 12, 14, shown in FIG. 1, is maintained. The width and thethickness of the stay 16 are respectively smaller than the length L_(W)and L_(H), of the first and second bars 12, 14. Further, a length alongthe longitudinal direction of stay 16 is shorter than the length of thelongitudinal direction of the first and second bars 12, 14.

The stay 16 is rotatably connected to the first bar 12 by a stay hinge92, being attachable to and removable from the second bar 14 by a lockhinge 94. When the target 10 is utilized in an operational positionshown in FIG. 1, the stay 16 is positioned so as to make a predeterminedangle with each of the first and second bars 12, 14, so that the firstand second bars 12, 14 make a right angle.

On the top surface of the target 10, namely, on the top surfaces of thebars 12, 14 and the controlling unit box 20, three standard pointmembers 31, 34 and 36, and three assistant point members 32, 33 and 35are mounted. The standard point member 31 (first standard point member)is a circular plate, the diameter of which is smaller than the widthL_(W) of the first and second bars 12, 14. The standard point members 34(second standard point member), 36 (third standard point member) and theassistant point members 32, 33 and 35 are identical to the standardpoint member 31.

The standard point member 31 and the assistant point members 32, 33 aremounted on the top surface 12 e of the first bar 12. The standard pointmember 34 is mounted on the top surface 20 e of the controlling unit box20. The assistant point member 35 and the standard point member 36 aremounted on the top surface 14 e of the second bar 14. The assistantpoint members 32 and 33 are positioned in such a manner that centers ofthe assistant point members 32 and 33 lie on a straight line (firststraight line) parallel to the axis 12 x, which connects a center of thestandard point member 31 and a center of the standard point member 34.Also, the assistant point member 35 is positioned in such a manner thata center of the assistant point member 35 lies on a straight line(second straight line) parallel to the axis 14 x, which connects acenter of the standard point member 34 and a center of the standardpoint member 36.

As described above, there are two assistant point members (32, 33) onthe straight line connecting the standard point members 31 and 34, andthere is one assistant point member (35) on the straight line connectingthe standard point members 34 and 36. Namely, the number of assistantpoint member on each of the straight lines is different from each other.Further, the distance between the center of the standard point member 31and the center of the assistant point member 32, the distance betweenthe center of the assistant point member 32 and the center of theassistant point member 33, and the distance between the center of theassistant point member 33 and the center of the standard point member34, are identical to each other. The distance between the center of thestandard point member 34 and the center of the assistant point member35, and the distance between the center of the assistant point member 35and the center of the standard point member 36, are identical to eachother. Furthermore, the distance between the center of the standardpoint member 31 and the standard point member 34, and the distancebetween the center of the standard point member 34 and the center of thestandard point member 36, are identical to each other.

The aforementioned reference plane for the photogrammetric analyticalmeasurement is defined by the standard point members 31, 34 and 36, andthe assistant point members 32, 33 and 35. The standard measurementlength is defined by a side length of an isosceles triangle apexes ofwhich correspond to the standard point members 31, 34 and 36. Namely, alength of the distance between the standard point members 34 and 36, anda length of the distance between the standard point members 36 and 31can be used the standard measurement length, as they are predetermined.

Note that, the angle θ is not restricted to 90°, and it is unnecessaryto make the distance between 31, 34 and the distance between 34, 36equal to each other. Namely, it is required that each value of the angleθ, the length of the distance between the standard point members 31, 34and the length of the distance between the standard point members 34, 36are predetermined. Considering facility of calculation in thephotogrammetric analytical measurement system, it is preferable that thevalue of the angle θ is predetermined to be 90° and the distance betweenthe standard point members 31, 34 and the distance between the standardpoint members 34, 36 are equal to each other.

As is apparent from FIG. 1, on each of the equilaterals of the isoscelestriangle, a number of the assistant point members is different.Accordingly, as a positional relationship between the target 10 andother objects in a picture can be easily judged, determination of thephotographing position of each picture for the photogrammetricanalytical measurement system is facilitated.

It is desirable that the standard plane defined by the standard pointmembers 31, 34, 36 and the assistant point members 32, 33, 35 isparallel to the surface of a road on which the target 10 is positioned.As the controlling unit box 20 is mounted on the first bar 12, the firstbar 12 is heavier than the second bar 14. Accordingly, if the second bar14 is connected to the first bar 12 only by the hinge 15 at the endportion thereof, the second bar 14 rises from the surface of the road,so that the standard point member 31 and the assistant point members 32,33 on the first bar 12 do not lie on the same plane on which theassistant point member 35 and the standard point member 36 on the secondbar 14 lie. Also, the standard plane might not be parallel to thesurface of the road.

Accordingly, the first and second bars 12, 14 are connected by the stay16 at center portions thereof, in order to prevent the second bar 14from rising from the surface of the road and to keep the standard planeparallel to the surface of the road for more accurate photogrammetricanalytic measurement.

The mounting of the hinge 15 produces an interval between the sidesurface 20 c of the controlling unit box 20 and the end surface 14 a ofthe second bar 14. An elastic member 19 (see FIG. 5) is disposed in theinterval. Accordingly, instability of the second bar 14 is prevented.The elastic member 19 is made of, for example, rubber and sponge, beingattached on the end surface 14 a or the side surface 20 c. Note that, aspring can be utilized instead of the elastic member 19.

A reflecting sheet is attached on the standard point members 31, 34, 36and the assistant point members 32, 33, 35. As the surface of thereflecting sheet is smooth, and colored white, the amount of reflectionof incident light on the sheet is high. The standard point members 31,34 and 36 are respectively surrounded by non-reflecting members 41, 44and 46. Similarly, the assistant point members 32, 33, 35 arerespectively surrounded by non-reflecting members 42, 43 and 45. On thenon-reflecting members 41, 42, 43, 44, 45 and 46, the non-reflectingsheet is attached, respectively. Accordingly, discrimination of thestandard point members 31, 34, 36 and the assistant point members 32,33, 35 in the pictures is facilitated, so that the photogrammetricanalytic measurement can be performed accurately.

The target 10 comprises a first tilt sensor 52 and a second tilt sensor54. Tilt angles of the target 10 around the axes 12 x and 14 x aresensed by the first and second tilt sensors 52, 54. The first tiltsensor 52 is mounted in the first bar 12, being placed between theassistant point members 32 and the standard point members 34. The tiltangle to a horizontal plane around the axis 12 x is sensed by the firsttilt sensor 52. The second tilt sensor 54 is mounted in the second bar14, being placed between the standard point member 34 and the assistantpoint member 35. The tilt angle of the horizontal plane around the axis14 x is sensed by the second tilt sensor 54.

The first and second tilt sensors 52, 54 are connected to thecontrolling box 20 by a cable 17. Data of the tilt angles is transmittedto the controlling box 20 by the cable 17.

The angle of the standard plane to the horizontal plane is obtained bysensing the tilt angles around the axes 12 x and 14 x which areperpendicular to each other.

The target 10 comprises three legs 18. The legs 18 are mounted on thebottom surface opposite to the top surface on which the standard pointmembers 31, 34, 36 and the assistant point members 32, 33, 35 aremounted. Each of the legs 18 respectively correspond to the standardpoint members 31, 34 and 36. In a spot for photogrammetric analyticmeasurement, the target 10 is supported above a road with the intervalcorresponding to the height of the leg 18 therebetween. Accordingly, thetarget 10 is positioned so as to be parallel to general inclination ofthe road, without being effected by roughness of the surface of theroad.

With reference to FIGS. 3 and 4, the constructions of the assistantpoint member 35 and the non-reflecting member 45 are described below.FIG. 3 is a sectional view of the target 10, taken in the direction ofthe arrows substantially along the line IX—IX of FIG. 1. FIG. 4 is aplan view of the non-reflecting member 45, shown from the side of thesecond bar 14.

A magnet holder 62 is mounted on the top surface 14 e of the second bar14. A magnet 64, which is ring-shaped, is held in the magnet holder 62.The outer diameter of the magnet holder 62 approximately equals thewidth L_(W) of the second bar 14. The magnet 64 and the magnet holder 62are integratedly fixed to the second bar 14 by a screw 66. A reflectingsheet 68 is attached to a head 67 of the screw 66. The assistant pointmember 35 comprises the magnet holder 62, the magnet 64, the screw 66and the reflecting sheet 68.

The non-reflecting member 45 comprises a circular plate 72, made of amaterial through which an electric (radio) wave transmits, for example,resin or rubber material. If the circular plate 72 is made of rubbermaterial, breakage of the circular plate 72 is avoided when the circularplate 72 is accidentally dropped. A non-reflecting sheet 74 is attachedto one surface of the circular plate 72. The diameter of thenon-reflecting member 45 is approximately seven times the diameter ofthe head 67 of the screw 66. The thickness of the non-reflecting member45 is slightly smaller than the thickness of the head 67.

An engagement hole 76 is formed at the center of the non-reflectingmember 45. The diameter of the engagement hole 76 approximately equalsthe diameter of the head 67. An iron ring 78 abuts around the engagementhole 76, at the surface opposite to the surface on which thenon-reflecting sheet 74 is attached. The inner diameter of the iron ring78 approximately equals the diameter of the engagement hole 76, and theouter diameter of the iron ring 78 approximately equals the outerdiameter of the magnet holder 62.

The non-reflecting member 45 is attachable to and removable from theassistant point member 35. When the target 10 is utilized for thephotogrammetric analytic measurement, the non-reflecting member 45 ismounted on the assistant point member 35 in such a manner that the head67 of the screw 66 is engaged with the engagement hole 76. In accordancewith the engagement of the head 67 and the engagement hole 76, the ironring 78 is fixedly attached to the magnet holder 62 by magnetic force ofthe magnet 64, so that the iron ring 78 is securely engaged with thehead 67. As is apparent from FIG. 3, when the non-reflecting member 45is mounted on the assistant point member 35, the reflecting sheet 68 andthe non-reflecting sheet 74 lie on a same plane. On the other hand, whenthe target 10 is not utilized, the non-reflecting member 45 is removedfrom the assistant point member 35 by a user, by resisting the magneticattraction between the magnet 64 and the iron ring 78.

Note that, the standard point members 31, 34, 36 and the assistant pointmembers 32, 33 are identical to the assistant point member 35. Further,the non-reflecting members 41, 42, 43, 44 and 46 are identical to thenon-reflecting member 45. Namely, each of the non-reflecting members isattachable to or removal from the corresponding standard point member.As each of the non-reflecting members is removable from thecorresponding standard point member, the target 10 is handy to carry.

Further, when the target 10 is utilized, the non-reflecting members arerespectively mounted on the corresponding standard point members, sothat each of the reflecting sheets (68) is surrounded by thenon-reflecting sheet (74). Photogrammetric analytic measurement may beperformed under a situation in which an amount of luminance is not idealfor photographing, for example, in rain or at nighttime. Also, themeasurement may be performed at a spot, in which a surface of a roadreflects incident light. However, even if photographing is performedunder such adverse conditions, areas of the standard point members 31,34, 36 and the assistant point members 32, 33, 35 can be viewed clearlyin photographed pictures, as each of the reflecting sheets is surroundedby the non-reflecting sheet.

Note that, the ratio between the diameter of the assistant point member35 and the diameter of the non-reflecting member 45, namely an amount ofarea of each of the reflecting sheet 68 and the non-reflecting sheet 74,is not restricted to this embodiment. It is just required that thereflecting sheet 68 is large enough for being able to be viewed clearlyin photographed pictures. Further, the assistant point member 35 and thenon-reflecting member 45 may have any shape other than a circle.

FIG. 5 shows an enlarged view of the controlling unit box 20 and anyother members close to the box 20, with portions broken away forclarity. FIG. 6 is a sectional view of the controlling unit box 20,taken in the direction of the arrows substantially along the lineXII—XII of FIG. 5. In FIG. 6, the construction of the controlling unitbox 20 is simplified.

The controlling unit box 20 comprises a battery room 83. The batteryroom 83 is positioned at the side of the side surface 20 d. A battery87, which supplies an electric power to the target 10, is installed inthe battery room 83. The battery room 83 has an opening in the side ofthe side surface 20 d. The opening is closed by a cover 83 a. A switch85 is unitarily mounted on the side surface 20 d. The power supply tothe target 10 is stopped and started by manipulation of the switch 85.

An opening 81 is formed in the top surface 20 e of the controlling unitbox 20. The opening 81 is closed by a cover 82. The cover 82 is made ofa material through which an electric (radio) wave can transmit, forexample, resin material. An antenna 88 is mounted on an inner surface ofthe cover 82, in such a manner that the antenna 88 is wound along theperiphery of the cover 82. In the controlling unit box 20, a controlboard 84 and an azimuth sensor 86 are mounted. The azimuth sensor 86 andthe tilt sensors 52, 54 are connected to the control board 84, wherebythe operation of the sensors is controlled by the control board 84.

As is apparent from FIG. 5, the azimuth sensor 86 is placed between theassistant point member 33 and the standard point member 34. Namely, whenthe distance between the assistant point member 33 and the standardpoint member 34 is L_(A), the distance between the center of the azimuthsensor 86 and the assistant point member 33 is L_(A)/2.

The sensing of the azimuth sensor 86 is affected by a magnetic materialplaced close to the sensor 86, for example, the controlling unit box 20and the battery 87. Further, as the non-reflecting members 41, 42, 43,44, 45 and 46 are respectively attached to the members 31, 32, 33, 34,35 and 36 by magnetic force, by which the sensing of the azimuth sensor86 may be affected. If the sensing of the azimuth sensor 86 is affectedby the magnetic force, it is necessary to compensate azimuth data sensedby the azimuth sensor 86. However, as described above, the azimuthsensor 86 is intermediately placed between the assistant point member 33and the standard point member 34. Namely, at the portion where thesensor 86 is placed, an influence of the magnetic force generated aroundthe assistant point member 33 and the standard point member 34 is at aminimum. Accordingly, the sensing of the azimuth sensor 86 is notsignificantly affected.

When the switch 85 is turned ON to start the power supply, the sensors52 and 54 measure the tilt angles and the azimuth sensor measures theazimuth at regular intervals in accordance with a control pulse outputfrom the control board 84. Data measured by the sensors 52, 54 and 86 isoutput to the control board 84. After the data is subjected topredetermined operations, for example, compensation, in the controlboard 84, the data is transmitted to a receiver from the antenna 86 bywireless. The receiver is mounted on an external device, for example, adigital camera (not shown).

As described above, the antenna 88 is placed on the inner surface of thecover 82. As the cover 82, the non-reflecting member 43, 44 are made ofthe material through which an electric (radio) wave can transmit, theelectric wave output from the antenna 88 is transmitted to the receiver,without being interrupted by the cover 82 and the non-reflecting member43, 44.

As described above, if the receiver is mounted in a digital camera, thetilt angle data, the azimuth and image data photographed by the digitalcamera are saved to a recording medium. By inputting such data to acomputer from the recording medium, processing of the image data can beperformed in the computer more rapidly, so that a survey map can be mademore accurately.

A spot of the photogrammetric analytic measurement, with the target 10placed at a predetermined position on the spot, is photographed by acamera,.

A procedure for calculating positional data of a camera using the target10 in this embodiment will now be explained. FIG. 7 shows a positionalrelationship between the target 10 and a camera 100. The target 10 isphotographed by the camera 100 at a camera position M. The cameraposition M is defined as a back principal point of a photographing lenssystem of the camera 100. An optical axis of the photographing lenssystem is represented by a broken line O. As described above, the target10 is provided with the three standard point members and the threeassistant point members thereon. However, only the three standard pointmembers are used for clarity of explanation. In FIG. 7, a plane, definedby the standard point members 31, 34 and 36, is the standard plane. Thelength of the distance between the standard point members 31, 34, andthe length of the distance between the standard point members 34, 36 arethe standard length L.

Note that, as is apparent from FIG. 7, an IC card 24, a readablerecording medium, is mounted in the camera 100. The IC card 24 isremovable from the camera 100, and attachable to other devices. Imagedata of the photographed pictures is recorded in the IC card 24. Whenthe photogrammetric analytic measurement is performed, the image datarecorded in the IC card 24 is loaded in a memory of the device of thephotogrammetric analytic measurement as described below, after mountingthe IC card 24 in the device of the photogrammetric analyticmeasurement.

FIG. 8 shows a photographed image at the camera position M. Atwo-dimensional coordinate system (Xp, Yp) is set in the photographedimage. An origin of the two-dimensional coordinate system is aphotographing center point CP. As apparent from FIG. 2, in thephotographed image, image points of the standard point members 31, 34and 36 respectively correspond to two-dimensional coordinates P1(Xp₁,Yp₁), P2(Xp₂, Yp₂) and P3(Xp₃, Yp₃). Note that, the two-dimensionalcoordinate system of the photographed image is referred to as a“photograph coordinate system”.

FIG. 9 shows a positional relationship between the photographed imagephotographed by the camera 100 and the target 10. A hatched area is thestandard plane defined by the standard point members 31, 34 and 36. InFIG. 9, a three-dimensional coordinate system of the camera position Mis shown as (Xc, Yc, Zc). An origin of the three-dimensional coordinatesystem corresponds to the back principal point of the photographing lenssystem of the camera 100. Further, the Zc-axis corresponds to theoptical axis O of the photographing lens system at the camera positionM. The Xc-axis is parallel to the Xp axis of the two-dimensionalcoordinate system of the photographed image shown in FIG. 8, and theYc-axis is parallel to the Yp axis of the two-dimensional coordinatesystem. In this specification, the three-dimensional coordinate systemat the camera position M is referred to as a “camera coordinate system”.

A three-dimensional coordinate system (Xs, Ys, Zs) shown in FIG. 10,which is a right-handed coordinate system, is set in order to determinethe camera position M based on the photographed image. An origin of thethree-dimensional coordinate system (Xs, Ys, Zs) is accorded to thestandard point member 34 of the target 10. A Zs-axis is a direction fromthe standard point member 34 to the standard point member 36. An Xs-axisis perpendicular to the Zs-axis, lying on the standard plane. A Ys-axisis an axis on which the standard point member 34 lies, beingperpendicular to the standard plane, namely being perpendicular to thesheet of drawing. As described above, the length of the distance betweenthe standard point members 31, 34 and the length of the distance betweenthe standard point members 34, 36 are L. Accordingly, image points ofthe standard point members 31, 34 and 36 respectively correspond tothree-dimensional coordinates Ps1(−L, 0, 0), Ps2(0, 0, 0) and Ps3(0, 0,L). In this specification, the three-dimensional coordinate system ofFIG. 10 is referred to as a “scene coordinate system”.

The camera position M is defined by a relationship between the scenecoordinate system and the camera coordinate system. Namely, the cameraposition M is defined by a relationship of a movement distance of theorigin of the camera coordinate system from the origin of the scenecoordinate system to, a rotational angle of the Xs-axis around theXc-axis, a rotational angle of the Ys-axis around the Yc-axis and arotational angle of the Zs-axis around the Zc-axis.

A relationship between Psi(Psxi, Psyi, Pszi) (i=1, 2, 3) in the scenecoordinate system and Pci(Pcxi, Pcyi, Pczi) (i=1, 2, 3) in the cameracoordinate system is stated in an equation (1) as below.

$\begin{matrix}{{P\quad c\quad i} = {{R\left( {{P\quad s\quad i} - \Delta} \right)}:}} & (1) \\{{{P\quad c\quad i} = \begin{pmatrix}{P\quad c\quad x\quad i} \\{P\quad c\quad y\quad i} \\{P\quad c\quad z\quad i}\end{pmatrix}};\quad {{P\quad s\quad i} = \begin{pmatrix}{P\quad s\quad x\quad i} \\{P\quad s\quad y\quad i} \\{P\quad s\quad z\quad i}\end{pmatrix}};} & \quad \\{{R = \begin{pmatrix}{{{Cos}\quad \beta \quad {Cos}\quad \gamma}\quad} & {{{Cos}\quad \alpha \quad {Sin}\quad \gamma} + {{Sin}\quad \alpha \quad {Sin}\quad \beta \quad {Cos}\quad \gamma}} & {{{{Sin}\quad \alpha \quad {Sin}\quad \gamma}\quad - {{Cos}\quad \alpha \quad {Sin}\quad \beta \quad {Cos}\quad \gamma}}\quad} \\{{{- {Cos}}\quad \beta \quad {Cos}\quad \gamma}\quad} & {{{- {Cos}}\quad \alpha \quad {Cos}\quad \gamma} - {{Sin}\quad \alpha \quad {Sin}\quad \beta \quad {Sin}\quad \gamma}} & {{{Sin}\quad \alpha \quad {Cos}\quad \gamma} + {{Cos}\quad \alpha \quad {Sin}\quad \beta \quad {Sin}\quad \gamma}} \\{{Sin}\quad \beta} & {{- {Sin}}\quad \alpha \quad {Cos}\quad \beta} & {{{Cos}\quad \alpha \quad {Cos}\quad \beta}\quad}\end{pmatrix}};} & \quad \\{{\Delta = \begin{pmatrix}{\Delta \quad X} \\{\Delta \quad Y} \\{\Delta \quad Z}\end{pmatrix}}\left( {{i = 1},2,3} \right)} & \quad\end{matrix}$

In the equation (1), α is the rotational angle of the Xs-axis around theXc-axis, β is the rotational angle of the Ys-axis around the Yc-axis, γis the rotational angle of the Zs-axis around the Zc-axis. Further, R isa rotational matrix of α, β and γ, and Δ is a vector from the scenecoordinate system to the camera coordinate system.

A relationship between Pci (Pcxi, Pcyi, Pczi) and Di (Dxi, Dyi) (i=1, 2,3) is stated in an equation (2) as below. Note that, Pci is thethree-dimensional coordinates of the standard point Pi in the cameracoordinate system, and Di is the two-dimensional coordinates of thestandard point in the photograph coordinate system. Futher, in theequation (2), “f” is a focal distance of the photographing lens systemof the camera 10.

Dxi=f·Pcxi/Pczi Dyi=f·Pcyi/Pczi}  (2)

(i=1,2,3)

As described above, the two-dimensional coordinates of each of thestandard points in the photograph coordinate are calculated from thethree-dimensional coordinates in the scene coordinate system based onthe equations (1) and (2). On the other hand, the two-dimensionalcoordinates of each of the standard points are automatically extractedfrom the photographed image. Accordingly, with respect to the standardpoints, by comparing the values of the two-dimensional coordinatescalculated from the equations (1) and (2) with the values of thetwo-dimensional coordinates extracted from the photographed image, thevalues defining the camera position, namely ΔX, ΔY, ΔZ, α, β and γ canbe calculated.

FIG. 17 shows a block diagram of a coordinate-calculator system, inwhich the photogrammetric analytic measurement, as mentioned above, isperformed on the basis of the digital image data stored in the IC card24.

As shown in FIG. 17, the coordinate-calculator system is constructed asa computer system comprising a computer 184, which includes: a centralprocessing unit (CPU) 184A; a read-only-memory (ROM) 184B having anoperating system program, constants, etc. stored therein; arandom-access-memory (RAM) 184C for storing temporary data, temporaryconstants, etc.; a hard disk 184D for storing data resulting fromcalculations executed by the CPU 184A; and an input/output interface(I/O) 184E.

The computer system also comprises an IC memory card reader 186connected to the computer 184, via the I/O 184E. The IC memory-cardreader 186 is provided with a slot for receiving the IC memory card 24,and includes an IC card driver 188 for reading a given frame of digitalimage data and other information data.

The computer system also comprises a CD-ROM reader 197 connected to thecomputer 184, via the I/O 184E. The CD-ROM reader 197 is provided with atray for receiving a CD-ROM 198, and includes a CD-ROM driver unit 199.A calculating program of the standard points of the target andpositional data of the camera, as described below, is stored in theCD-ROM 198. The calculating program is read from the CD-ROM 198, andinstalled to the hard disk 184D.

The computer system further comprises a monitor 190 for reproducing aphotographed picture based on the frame of digital image data read fromthe IC memory card 24 and a survey map produced by the computer 184, akeyboard 192 for inputting various command signals and various data tothe computer 184, a mouse 194 for manipulating a cursor displayed on themonitor 190, and a printer 196 for printing the survey map on a sheet ofpaper, if necessary.

FIG. 11 shows a flowchart indicating a process of the calculatingprogram of the camera position, according to this embodiment, executedby the CPU 184A.

After the target 10 is situated on a predetermined position in a spotfor photogrammetric analytic measurement, the spot is photographed bythe camera 100. In step S1000, digital image data of the spot is loadedto a system controller of a photogrammetric analytic measurement deviceto which this embodiment is applied.

In step S1002, a gray scale operation is performed. By the gray scaleoperation, brightness information is calculated based on colorinformation of one frame color image data. After the brightnessinformation, which is represented by 256 intensity levels, with respectto each pixel composing the digital image data, the brightnessinformation is stored to a variable Img of array type. Namely, thebrightness information of each pixel is represented by Img (j, k). Notethat, “j ” and “k” correspond to the values of the two-dimensionalcoordinates (the photograph coordinate system) of each pixel in thedigital image data.

In step S1004, binarization of the brightness information of each pixelis performed. In this embodiment, the binarization is performed usingthresholding. The brightness information stored in the variable Img iscompared with a predetermined threshold. If the brightness informationis higher than the predetermined threshold, “1” is stored to a variableBin of array type. If the brightness information is not higher than thepredetermined threshold, “0” is stored to the variable Bin. Namely,binary information of the brightness information of each pixel isrepresented by Bin (j, k).

As described above, the reflecting sheet is attached on the standardpoint members 31, 34, 36 and the assistant point members 32, 33, 35. Thestandard point members 31, 34 and 36 are respectively surrounded by thenon-reflecting members 41, 44 and 46 which are black. Likewise, theassistant point members 32, 33, 35 are respectively surrounded by thenon-reflecting members 42, 43 and 45 which are black. Accordingly, thebinary information of pixels corresponding to the standard point members31, 34, 36 and the assistant point members 32, 33, 35 is represented by“1”, and the binary information of pixels corresponding to thenon-reflecting members 41, 42, 43, 44, 45, 46 is represented by “0”.

Then, in step S1006, labeling is performed. By the labeling, a set ofpixels, binary information of which are “1”, is extracted as a group.Namely, regions, brightness information of which are high, are extractedfrom the digital image data.

In step S1008, an area of each group is calculated. Namely, with respectto each group, a number of pixels composing the group is calculated. Thenumber of pixels of each group and two-dimensional coordinatescorresponding to the pixels are stored to predetermined variables.

Then, in step S1010, it is judged based on the number of pixels of eachgroup if each group corresponds to the standard point members or theassistant point members. It is checked that the number of the pixels isbetween a greatest lower bound Th1 and a least upper bound Th2. Thevalues of Th1 and Th2 are determined, such that a range of the pixelnumbers defined by Th1 and Th2 is able to contain the pixel number ofthe group corresponding to the standard point members and the assistantpoint members, with the range being as minimal as possible.

Accordingly, a group corresponding to noise of the digital image, apixel number of which is one or two, and a group corresponding to a wideregion of high brightness, for example, a painted region of bright coloron a road of the spot, are deleted.

In step S1012, a center of gravity (spectral center) of each group iscalculated based on the brightness information of each group. FIG. 12shows a flowchart of a routine for calculating the center of gravity. Instep S1020, one group is selected. In step S1022, variables SumX, SumYand SumImg are initialized. A continued product of the brightnessinformation of each pixel composing the selected group and an Xcomponent of the two-dimensional coordinate of the each pixel is storedto the SumX. A continued product of the brightness information of eachpixel composing the group and a Y component of the two-dimensionalcoordinate of the each pixel is stored to the SumY. A sum total of thebrightness information of each pixel composing the group is stored tothe SumImg.

In step S1026, with respect to one pixel composing the group, a productof the brightness information (Img (j, k)) and the X component (j) iscalculated to be stored to the SumX after being added to a current valuestored in the SumX. Then, in step S1028, with respect to the same pixelcomposing the group, a product of the brightness information (Img (j,k)) and the Y component (k) is calculated to be stored to the SumY afterbeing added to a current value stored in the SumY. In step S1030, thebrightness information of the same pixel is added to a current valuestored in the SumImg to be stored to the SumImg. The above mentionedprocesses from step S1026 through step S1030 are performed repeatedly,with respect to each pixel composing the group.

Namely, a sum total of weighted brightness information based on Xcomponent of each pixel composing the group is stored in the SumX, and asum total of weighted brightness information based on Y component ofeach pixel composing the group is stored in the SumY. As describedabove, the sum total of the brightness information of the pixels of thegroup is stored in the SumImg.

In step S1032, it is checked if all pixels composing the group have beensubjected to the above mentioned procedures. If the above mentionedprocedures have been performed with respect to all pixels composing thegroup, the process goes to step S1034. In step S1034, an X component ofthe center of gravity is calculated by dividing the sum total ofweighted brightness information based on X component of each pixel bythe sum total of the brightness information of each pixel. In stepS1036, a Y component of the center of gravity is calculated by dividingthe sum total of weighted brightness information based on Y component ofeach pixel by the sum total of the brightness information of each pixel.

Then, in step S1038, it is checked if all groups in the digital imagedata have been subjected to the calculation of the center of gravity. Ifthere are any groups which have not been subjected to the calculation,the procedures from step S1020 through step S1036 are repeatedlyperformed. If all groups have been subjected to the calculation of thecenter of gravity, the process goes to step S1014 of FIG. 11.

FIG. 13 is a conceptual view of each pixel of the digital image dataafter the calculation of the center of gravity. In FIG. 13, pixels oflow brightness information are hatched by oblique lines. As is apparentfrom FIG. 13, groups from G1 through G13 have been extracted. Thecenters of gravity of the groups are respectively represented by GC1through GC13.

In step S1014, combinations of the groups, which satisfy a positionalrelationship between the standard position members 31, 34, 36 and theassistant position members 32, 33, 35 of the target 10, are extracted.FIG. 14 shows the first half of a flowchart of an extraction routine ofthe combinations of the groups. FIG. 15 shows the latter half of theflowchart. The procedures of the flowchart will be explained withreference to FIG. 13.

In step S1040, two groups are optionally selected. In step S1042, astraight line, which lies on the centers of gravity of the selected twogroups, is calculated. Then, in step S1044, groups which have a centerof gravity between the centers of gravity of the selected two groups,lying on the straight line, are detected.

In step S1046, a number of the groups detected in step S1044 is checked.If there are four groups, which have the center of gravity on thestraight line, including the selected two groups in step S1040, theprocess goes to step S1048. If there are three groups, which have thecenter of gravity on the straight line, including the selected twogroups in step S1040, the process goes to step S1050. Otherwise, theprocess goes to step S1052.

For example, if the group G1 and the group G2 are selected in stepS1040, there is no other group which has a center of gravity between thecenters of gravity GC1 and GC2 on a straight line on which the center ofgravity GC1 and the center of gravity GC2 lie. Accordingly, the processgoes to step S1052. If the groups G3 and G6 are selected, there are twocenters of gravity, GC4 of a group G4 and GC5 of a group G5, on astraight line on which a center of gravity GC3 of the group G3 and acenter of gravity GC6 of the group G6 lie. Accordingly, the process goesto step S1048.

Also, if the groups G10 and G13 are selected, there are two centers ofgravity, GC11 of a group G11 and GC12 of a group G12, on a straight lineon which a center of gravity GC10 of the group G10 and a center ofgravity GC13 of the group G13 lie. Accordingly, the process goes to stepS1048.

If the groups G1 and G3 are selected, there is one center of gravity,GC2 of the group G2, on a straight line on which the center of gravityGC1 of the group G1 and a center of gravity GC3 of the group G3 lie.Accordingly, the process goes to step S1050.

If the groups G7 and G9 are selected, there is one center of gravity,GC8 of a group G8, on a straight line on which a center of gravity GC7of the group G7 line. Accordingly, the process goes to step S1050.

In step S1048, the four groups are registered as a four points rangewhich lie on a same straight line. In step S1050, the three groups areregistered as a three points range which lie on a same straight line.For example, a set of the groups G1, G2, G3 and a set of the groups G6,G7, G8 are respectively registered as three points ranges. Also, a setof the groups G4, G5, G6, G7 and a set of the groups G10, G11, G12, G13are respectively registered as four points ranges.

In step S1052, it is checked if all combinations of two groups in thedigital image data have been subjected to the procedures from step S1042through step S1050. If there are some combinations of two groups whichhave not been subjected to the procedures, the process returns to stepS1040 and the procedures from step S1040 through step S1050 arerepeatedly performed. If there is no combination of two groups whichhave not been subjected to the procedures, the process goes to stepS1054 of a flowchart of FIG. 15. This process may be considered to besomewhat flexible. That is, the process for detection of gravity on astraight line between the centers of gravity of the two selected groupsshould be considered to detect centers of gravity substantially on thestraight line.

In step S1054, one combination of the registered three points range andthe registered four points range is selected. Then, in step S1056, it ischecked if a center of gravity of a group corresponding to end point ofthe registered three points range and a center of gravity of a groupcorresponding to end point of the registered four points range are inaccord, with respect to the combination selected instep S1054. Namely,it is checked if there is an identical end point shared by theregistered three points range and the registered four points range, ofthe combination. If the end points accord, the process goes to stepS1058. If the end points do not accord, namely, the registered threepoints range and the registered four points do not share an identicalend point, the process returns to step S1054 to select anothercombination of the registered three points range and the registered fourpoints range. In FIG. 13, the group G3 is an endpoint of the registeredthree points range composed by the groups G1, G2, G3, and also an endpoint of the registered four points range composed by the groups G3, G4,GS, G6. Accordingly, if the registered three points range of the groupsG1, G2, G3 and the registered four points range of the groups G3, G4,G5, G6 are selected in step S1054, the process goes to step S1058.

In step S1058, two-dimensional coordinate of a center of gravity of thegroup corresponding to another end point, of the registered four pointsrange, opposite to the end point which accords with the end point of theregistered three points range, is stored to a variable A1 (Ax₁, Ay₁).Also, two-dimensional coordinate of the center of gravity of the groupcorresponding to the end point shared by the registered three pointsrange and the registered four points range is stored to a variable A2(Ax₂, Ay₂). Further, two-dimensional coordinate of a center of gravityof the group corresponding to another end point, of the registered threepoints range, opposite to the end point which accords with the end pointof the registered four points range, is stored to a variable A3 (Ax₃,Ay₃). Namely, two-dimensional coordinates of the center of gravity GC6are stored to the variable A1, two-dimensional coordinates of the centerof gravity GC3 are stored to the variable A2 and, two-dimensionalcoordinates of the center of gravity GC1 are stored to the variable A3.

In step S1060, an angle θ, by which a vector from the variable A2 to thevariable A3 is rotated in the counterclockwise direction from a vectorfrom the variable A2 to the variable A1, is calculated. It is checkedthat the angle θ is larger than a greatest lower bound Th3, beingsmaller than a least upper bound Th4. Note that, in this embodiment, 0degrees is set to the greatest lower bound Th3, and 180 degrees is setto the least upper bound Th4. If the angle θ is between the greatestlower bound Th3 and the least upper bound Th4, the selected combinationof the registered three points range and the registered four pointsrange meets a positional relationship between the standard point members31, 34, 36 and the assistant point members 32, 33, 35. Accordingly, theselection of the combination is ended and the process goes to step S1016of the flowchart in FIG. 11.

If the angle θ is smaller than the greatest lower bound Th3 or largerthan the least upper bound Th4, the process goes to step S1062. In stepS1062, it is checked if all combinations of the registered three pointsrange and the registered four points range have been subjected to theprocedures from step S1056 through step S1060. While there are somecombinations which have not been subjected to the procedures from stepS1056 through step S1060, the procedures are repeatedly performed.

Namely, until a combination of the registered three points range and theregistered four points range, which meets the positional relationshipbetween the standard point members 31, 34, 36 and the assistant pointmembers 32, 33, 35, is found, the procedures from step S1054 throughstep S1062 are repeatedly performed. When one combination, which meetsthe positional relationship between the standard point members and theassistant point members, is found, the procedures from step S1054through step S1062 cease to be performed and the selection ofcombination of the registered three points range and the registered fourpoints range is ended.

As described above, with reference to FIG. 13, a set of the groups G1,G2, G3 is selected as the three points range, and a set of the groupsG3, G4, G5, G6 is selected as the four points range, and thetwo-dimensional coordinates of GC6 is stored to the variable A1, thetwo-dimensional coordinates of GC3 is stored to the variable A2 and, thetwo-dimensional coordinates of GC1 is stored to the variable A3. As isapparent from FIG. 13, an angle by which the vector from the variable A2to the variable A3 is rotated in the counterclockwise direction from thevector from the variable A2 to the variable A1, is between the greatestlower bound Th3 and the least upper bound Th4. Accordingly, thecombination of the groups G1, G2, G3 and the groups G3, G4, G5, G6 isextracted as a combination which meets the positional relationshipbetween the standard point members 31, 34, 36 and the assistant pointmembers 32, 33, 35.

In step S1016, a calculation of a camera position is performed. FIG. 16shows a flowchart of a calculation routine of the camera position. Instep S1070, variables ΔX, ΔY, ΔZ, α, β and γ of the aforementionedequation (1) are initialized.

In step S1072, (−L, 0, 0), (0, 0, 0), (0, 0, L) are respectivelysubstituted for Ps1, Ps2, Ps3 of the aforementioned equation (1), andsuitable values are respectively substituted for the variables ΔX, ΔY,ΔZ, α, β and γ. Then, three-dimensional coordinates of the standardpoint members 31, 34 and 36 in the camera coordinate system, namelyPc1(Pcx₁, Pcy₁, Pcz₁), Pc2(Pcx₂, Pcy₂, Pcz₂) and Pc3(Pcx₃, PCy₃, Pcz₃)are calculated.

In step S1074, three-dimensional coordinates of Pc1, Pc2, Pc3 calculatedin step S1072, are substituted for Pc1, Pc2 and Pc3 of theaforementioned equation (2). Then, two-dimensional coordinates ofD1(Dx₁, Dy₁), D2(Dx₁, Dy₁), D3(Dx₁, Dy₁) in the photograph coordinatesystem are calculated.

In step S1076, two-dimensional coordinates of A1, A2, A3 in thephotograph coordinate system, which are automatically extracted from thedigital image data, are compared with D1, D2, D3 calculated from theequations (1) and (2). Namely, an absolute value of difference betweenD1 and A1, an absolute value of difference between D2 and A2, anabsolute value of difference between D3 and A3 are calculated, and thesum S of the absolute values is calculated.

In step S1078, it is checked if the sum S is at a minimum. If the sums Sis not at a minimum, the procedures from step S1072 through step S1076are repeated after each value of ΔX, ΔY, ΔZ, α, β, γ is changed. Namely,values of ΔX, ΔY, ΔZ, α, β, γ, by which the sum S is at a minimum, arecalculated by least squares method.

If it is judged that the sum S is at a minimum in step S1078, theprocess goes to step S1080. Instep S1080, the values of ΔX, ΔY, ΔZ, α,β, γ are registered as parameters which define the positional data ofthe camera, and the calculation of the positional data of the camera isended.

The above mentioned procedures of the flowcharts shown in FIGS. 11, 12,14, 15 and 16 is the calculation of the two-dimensional coordinates ofthe standard points of the target and the positional data of the camera,in the photogrammetric analytic measurement. Namely, the above mentionedprocedures correspond to, for example, procedures of flowcharts fromstep 1201 through step 1211 shown in FIG. 12 of U.S. application Ser.No. 09/017,237. Due to the above mentioned procedures, calculation ofthree-dimensional coordinates of reference points (correspond to thestandard points in this embodiment) of step 1211 can be accuratelyperformed without being affected by the skill of an operator, asdesignation of the reference points of step 1208 is automaticallyperformed.

As described above, according to this embodiment, the number of theassistant point members between the standard point members 31 and 34,and the number of the assistant point members between the standard pointmembers 34, 36, are different. Further, the reflecting sheet is attachedon the standard point members and the assistant point members.Accordingly, it is easy to determine the positions of the standard pointmembers and the assistant point members, in the photographed picture, sothat two-dimensional coordinates of the standard point members in thedigital image are automatically calculated.

Further, according to this embodiment, as the length between thestandard point members 31, 34, the length between the standard pointmembers 34, 36, and the angle defined by the first bar 12 and the secondbar 14 are determined, the three-dimensional coordinates of the standardpoint members in the scene coordinate system are predetermined. Further,the three-dimensional coordinates of the standard point members in thephotograph coordinate system are calculated based on thethree-dimensional coordinates in the scene coordinate system.Accordingly, the positional data of the camera is automaticallycalculated by comparing the two-dimensional coordinates of the standardpoint members in the photograph coordinate system with thetwo-dimensional coordinates of the standard point members in the digitalimage.

As described above, according to the present invention, the positionaldata of a camera is automatically calculated in the photogrammetricanalytic measurement.

The present disclosure relates to subject matter contained in thefollowing Japanese Patent Application No. 10-334320 (filed on Nov. 25,1998), which is expressly incorporated herein, by reference, in itsentirety.

What is claimed is:
 1. A device for calculating positional data ofstandard points of a target of photogrammetric analytic measurementcomprising: an image data obtaining processor that photographs an objectof photogrammetric analytic measurement with a target that has a firststandard point and a second standard point and a third standard point,each length of distance between said first standard point and saidsecond standard point and said third standard point being predetermined,an angle defined by a first straight line connecting said first standardpoint and said second standard point, and a second straight lineconnecting said second standard point and said third standard point,being predetermined, a number of assistant points provided on said firststraight line and a number of assistant points provided on said secondstraight line being a different number of assistant points; and a valueof two-dimensional coordinates calculating processor that identifies aposition of each of said first standard point, said second standardpoint and said third standard point in said photographed image inaccordance with said different number of assistant points provided onsaid first straight line and provided on said second straight line; andthen calculates each two-dimensional coordinate value of said firststandard point, said second standard point and said third standard pointin a photograph coordinate system that is a two-dimensional coordinatesystem of a photographed image, obtained by said image data obtainingprocessor.
 2. The device of claim 1, wherein said first standard point,said second standard point, said third standard point and said assistantpoints provided on said first straight line and provided on said secondline are made of a reflecting material, such that said first standardpoint, said second standard point, said third standard point and saidassistant points have a relatively higher brightness than other objectsin said photographed image, and said value of two-dimensionalcoordinates calculating processor is provided with: a binarizationprocessor that classifies each pixel composing of digital image data ofsaid photographed image into a high brightness class and a lowbrightness class after comparing a brightness of each pixel with apredetermined threshold value; a group extracting processor thatextracts groups in which a plurality of said pixels of said highbrightness class sequentially exist in said digital image data, an areaof which is within a predetermined range; a calculating processor thatcalculates a value of two-dimensional coordinates of a center of gravityof each of said groups, based on said brightness of each pixel composingsaid groups and a value of coordinates of said each pixel in saidphotograph coordinate system; a first points range extracting processorthat extracts a first points range on which a plurality of said centersof gravity lie, a number of said plurality of said centers of gravity onsaid first points range being equal to a total point number includingsaid first standard point, said second standard point and said number ofsaid assistant points on said first straight line; a second points rangeextracting processor that extracts a second points range on which aplurality of said centers of gravity lie, a number of said plurality ofsaid centers of gravity on said second points range being equal to atotal point number including said second standard point, said thirdstandard point and said number of said assistant points on said secondstraight line; a combination extracting processor that extracts acombination of said first points range and said second points range inwhich said first points range and said second points range share saidcenter of gravity of an identical end point; and a standard pointdetermining processor that determines a center of gravity of an endpoint of said first points range extracted by said first points rangeextracting processor, being opposite to said identical end point sharedby said first points range and said second points range, as said firststandard point, determines said center of gravity of said identical endpoint shared by said first points range and said second points range assaid second standard point, and determines a center of gravity of an endpoint of said second points range extracted by said extractingprocessor, being opposite to said identical end point shared by saidfirst points range and said second points range, as said third standardpoint.
 3. The device of claim 2, wherein said number of said assistantpoints provided on said first straight line is two, and said number ofsaid assistant points provided on said second straight line is one, sothat said first points range is a four points range in which fourcenters of gravity range on a same straight line in said photographedimage, and said second points range is a three points range in whichthree centers of gravity range on a same straight line in saidphotographed image.
 4. A method for calculating positional data ofstandard points of a target of photogrammetric analytic measurementcomprising: a first step that photographs an object of photogrammetricanalytic measurement with a target that has a first standard point and asecond standard point and a third standard point, by a camera, eachlength of distance between said first standard point and said secondstandard point and said third standard point being predetermined, anangle defined by a first straight line connecting said first standardpoint and said second standard point, and a second straight lineconnecting said second standard point and said third standard point,being predetermined, a number of assistant points provided on said firststraight line and a number of assistant points provided on said secondstraight line being a different number of assistant points; and a secondstep that identifies a position of said first standard point, saidsecond standard point and said third standard point in said photographedimage in accordance with said different number of assistant pointsprovided on said first straight line and provided on said secondstraight line, and then calculates each two-dimensional coordinate valueof said first standard point, said second standard point and said thirdstandard point in a photograph coordinate system that is atwo-dimensional coordinate system of a photographed image photographedby said camera.
 5. A recording medium in which a calculating program ofstandard points of a target of photogrammetric analytic measurement isstored, said target having a first standard point and a second standardpoint and a third standard point, each length of distance between saidfirst standard point and said second standard point and said thirdstandard point being predetermined, an angle defined by a first straightline connecting said first standard point and said second standardpoint, and a second straight line connecting said second standard pointand said third standard point, being predetermined, a number ofassistant points provided on said first straight line and a number ofassistant points provided on said second straight line being a differentnumber of assistant points, said calculating program comprising: avalues of two-dimensional coordinate calculating routine that identifiesa position of each of said first standard point, said second standardpoint and said third standard point in said photographed image inaccordance with said different number of assistant points provided onsaid first straight line and provided on said second straight line; andthen calculates each two-dimensional coordinate value of said firststandard point, said second standard point and said third standard pointin a photograph coordinate system that is a two-dimensional coordinatesystem of a photographed image, photographed in such a manner that anobject of photogrammetric analytic measurement is photographed togetherwith said target by a camera.